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I elaborate on my prediction that an indirect detection of cold dark matter (CDM) 
may be possible by observing the gravitational lensing effects of the CDM cusp caustics 
at cosmological distances. Cusps in the distribution of CDM are plentiful once density 
perturbations enter the nonlinear regime of structure formation. Caustic ring model of 
galactic halo formation provides a well defined density profile and a geometry near the 
cusps of the caustic rings. I calculate the gravitational lensing effects of the cusps in 
this model. As a pointlike background source passes behind a cusp of a cosmological 
foreground halo, the magnification in its image may be detected by present instruments. 
Depending on the strength of detected effect and the time scale of brightness change, it 
may even be possible to discriminate between the CDM candidates: axions and weakly 
interacting massive particles. 
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1. Introduction 

Caustics of light 1 are generically the enveloping surfaces of a family of light rays 
where the intensity is high. Such concentration of light can burn, hence these en- 
velopes are called "caustics." The caustic formation in the propagation of light is 
a common phenomenon in nature. Rainbows, scintillation of stars, shimmering of 
the sea and dancing patterns of bright lines that appear on the bottom of a swim- 
ming pool are among the most familiar examples. When the wavefront of sun light 
enters into water, it gets retarded proportional to the local thickness of the wavy 
surface. If the thickness is small, the amplitude of the wavefront is not affected 
much, hence, the density of light rays remains almost uniform just under the water 
surface. In a deeper level from the surface, however, the rays necessarily intersect 
each other, because the tiny deviations in the paths of the rays due to the refraction 
are amplified. They form cusp-shaped regions bounded by an enveloping surface of 
rays. Three light rays hit every point inside this region, whereas one light ray hits 



1 
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every point outside. The enveloping surfaces, where the number of light rays jumps 
by two, are the fold caustics of light. Cusps occur wherever the folds merge. The 
density diverges when a caustic is approached from the side with extra rays; see Sec- 
tion 2. Thus, a mild modulation in the wavefront spontaneously generates a strong 
enhancement in intensity some distance away, at caustic locations. Caustics move 
about but remain robust under stirrings of the water surface. 

Caustics are common in the propagation of light, because (i) photons are clas- 
sically collisionless, and (ii) flow of light from a distant (pointlike) source has zero 
velocity dispersion. Caustics of ordinary luminous matter is unusual, because or- 
dinary matter is normally collisionful. However, the fall of stars in and out of a 
galactic gravitational potential well is collisionless. Many elliptic galaxies are sur- 
rounded by ripples 2 in the distribution of light that are interpreted as caustics of 
luminous matter. 3 These ripples provide an existence proof of discrete flows and 
caustics in the infall of ordinary matter which constitutes only 4% of the content in 
the Universe. It is believed that 22% of the composition consists of nonluminous, 
collisionless particles with very small primordial velocity dispersion. 4 Such particles 
are called cold dark matter (CDM). The leading CDM candidates are axions and 
weakly interacting massive particles (WIMPs). Galaxies are surrounded by CDM 
that keeps falling into the gravitational potential wells of galaxies from all direc- 
tions and forms halos around their baryonic disks. Because the number densities 
of axions and WIMPs are huge, their infall onto isolated galaxies is regarded as a 
continuous flow. The number of flows, however, is discrete 5 ' 6 at any point in the 
halo and changes as a function of space and time. The enveloping surfaces where 
the number of flows jumps by two are the CDM caustics. Sharp discontinuities in 
the CDM density distribution occur at the locations of the caustics. Existence of 
caustics in such a CDM infall 7-11 can be taken for granted, once the existence of 
ripples 3 in the infall of ordinary mater is given. In a typical galaxy, any effect one 
can imagine that may erase CDM caustics, would have already erased these ripples. 
Indeed, there are several evidences for CDM caustics found in the halos of galaxies 
including the Milky Way. 12-14 

There are at least two sets of caustics in the halos of isolated galaxies: outer 
caustics and caustic rings. 8-10 Outer caustics are simple fold (A 2 ) catastrophes 
located on topological spheres enveloping the galaxy. Caustic rings, on the other 
hand, are closed tubes whose transverse cross-section is a closed line with three 
(dual) cusps, called a "tricusp." One of the cusps, the outer cusp, lies in the galactic 
plane, the other two are nonplanar; see Fig. 1. The triaxial tricusp is called elliptic 
umbilic (-D-4) catastrophe. 

To see the caustic formation 8-10 in the infall of CDM, consider the time evolution 
of CDM particles that are about to fall into the potential well of an isolated galaxy 
at time t\ for the first time in their history. They form a topological sphere with 
radius R\{t\) enclosing the galaxy called the "turnaround" sphere at time t\. After 
falling all the way through the disk of the galaxy, the particles of this turnaround 
sphere form a new sphere which reaches its maximum radius #2(^2) at time i 2 - 
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-^2(^2) < Ri(ti) because of the deepening of the potential well in the meantime. 
Note also that, due to the growth of the galaxy, the first turnaround radius Rxfo) at 




Fig. 1. Cross-section of a caustic ring in the case of axial and reflection symmetry, p and z 
are galactocentric cylindrical coordinates, p and q are the transverse dimensions in the p and z 
directions, respectively, a is the ring radius, p and q are taken of order a, for clarity. p,q -C a for 
actual caustic rings. 

time t2 is larger than Ri(t\). The sphere continues to oscillate in this way. Therefore, 
there are discrete number of flows in and out of the galaxy at any location in the 
halo. The local number of flows increases, as one approaches the galactic center 
from an arbitrary direction at a given time. First, in the outermost regions, there 
is one flow (particles falling in for the first time), then, at a certain location, the 
number jumps to three (particles falling in for the first time, particles falling out 
for the first time, particles falling in for the second time), and then to five, to seven, 
and so on, at other certain locations closer to the center. The boundary — which 
is a topological sphere — between the region with one (three, five, . . .) flow and 
the region with three (five, seven, . . .) flows is the first (second, third, . . .) outer 
caustic. Outer caustics are simple fold catastrophes. In the limit of zero velocity 
dispersion, the density diverges when the outer caustics are approached from the 
side with two extra flows; see Section 3. To see the formation of caustic rings with 
tricusp cross-section, on the other hand, a closer inspection to the time evolution 
of the infall sphere is needed. a Notice that, each time the turnaround sphere falls 
through the galaxy, it turns itself inside out. 8-10 The particles near the equator of 
the infall sphere carry too much angular momentum to reach the central part of the 
galaxy, and therefore, just before the sphere turns itself inside out, there exists a 
circle of points which are inside the sphere last. This circle is the locus of outer cusps 
that lie in the plane of the galaxy (see Fig. 1). Hence, caustic rings are located near 
where the particles with the most angular momentum in a given inflow reach their 



a See Ref. 10 for details. 
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distance of closest approach to the galactic center before going back out. 10 Because 
the CDM infall is continuous, the caustic ring in space-time associated with the 
infall sphere is a persistent feature in space. There is a caustic ring due to CDM 
particles falling through the galaxy for the first time, a caustic ring of smaller radius 
due to particles falling through for the second time, and so on. Like outer caustics, 
caustic rings are boundary surfaces where the number of flows jumps by two. Inside 
the caustic rings there are four flows, whereas there are two flows outside. 

Caustic structure in the dark matter infall may be revealed by observing their 
gravitational lensing effects. 15-19 The effects are largest near the cusps 19 which are 
common features in the halo CDM distribution. 13 In this paper, I review the lensing 
effects near the cusps of caustic rings. 19 The outline is as follows. In Section 2, 
caustics in the propagation of light are introduced. They provide a useful analog for 
understanding the properties of caustics in the infall of cold dark matter. The outer 
caustics and caustic rings of CDM are reviewed in Section 3. Their gravitational 
lensing effects are discussed in Section 4. Conclusions are summarized in Section 5. 

2. Caustics of Light 

When natural processes focus light, the caustics that are generated have a system- 
atic pattern. 1 The mathematical tool to investigate this systematics is the catas- 
trophe theory which analyzes the different ways that the critical points of gradi- 
ent mappings can merge to make higher singularities. In optics, the mappings are 
from an incident wave front to the space in which the wave is observed; i.e. three- 
dimensional (3D) space around the caustic or 2D space of a screen. The quantity 
that is extremized is the Fermat's potential and the critical points correspond to 
the light rays. The annihilation of two critical points at a location by coalescence 
with each other means loss of two rays contributing to a focus at that location. 

Caustics spontaneously arise when a forefront object distorts the propagation 
of light from a distant (pointlike) source. Although the light rays are parallel ini- 
tially, they necessarily cross each other after the distortion and produce caustics 
systemized as folds and cusps. If an observer places his/her eye inside the region 
bounded by two folds that merge at a cusp, light rays enter eye from three distinct 
directions, each corresponding to an image of the source. Every point outside the 
region, on the other hand, is hit by one ray only. A 2D cross-section of such a region 
is depicted in Fig. 2. Thus, as the point of observation crosses a fold caustic surface, 
the number of images jumps by two. If the eye moves from the inside across one of 
the folds, two of the images approach each other increasing in brightness, coalesce 
at the fold and then disappear outside. Meanwhile, the third image remains pretty 
much the same. If the eye moves towards the other fold, the third image coalesces 



b For example, gravity of massive objects like stars, focus the CDM flow producing a cone-shaped 
caustic surface — whose apex is a cusp — near the trajectories of maximum scattering angle. 20 
This effect can be generalized considering an isolated galaxy or even an isolated (super)cluster of 
galaxies as a massive object, as a whole. 
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Fig. 2. Distribution of light rays in the vicinity of a two dimensional cut of a region bounded 
by two folds merging at a cusp. The local curvature radius R of the fold at point P along the 
tangential direction is that of the circle depicted by dashed lines. 

with one of the other two images at the fold and then disappear. The closer the eye 
is to the cusp, the closer the two images are. The three images merge at the cusp. 
Because rays are more concentrated near the cusps, the intensity is largest there; 
see Fig. 2. When a fold is approached from the side with two extra rays at every 
point, the density rises as er~ 1//2 in the vicinity of the fold, where a is the distance to 
the fold in the direction of approach. This can be seen as follows. The rays tangent 
to a small segment of a smooth caustic fold at a point P are uniformly distributed 1 
along it; see Fig. 2. The corresponding length of the segment ~ 2\ / 2R<r where R is 
the local curvature radius of the fold at P along the tangential direction. Therefore, 
in the caustic neighborhood, the number of rays N passing between the observer 
and the fold is proportional to 2\/2Ro. Hence, the number density dN/da <~ cr -1 / 2 . 
Because caustics of light form spontaneously in rather general circumstances and 
retain their identity under perturbations, they are generic and structurally stable. 

3. Caustics of Cold Dark Matter 

CDM infall onto isolated galaxies intrinsically has the necessary conditions for the 
caustic formation noted in Section 1 : (i) the particles interact only weakly, and (ii) 
the primordial velocity dispersions Sv of the CDM candidates are of a very small 
order. For the axion infall, the primordial velocity dispersion at time t, is of order 8 ' 21 
5v a (t) <~ 10~ 12 km/s (lO~ 5 cV/r7i a ) 5 ^ 6 (to/t) 2 ^ 3 , where to an d m a are respectively 
the present age of the Universe and the axion mass. For the WIMP infall, Sv x (t) ~ 
10- 6 km/s (GeV/m x ) 1/2 (McV /T D ) 1/2 (t a /t) 2/3 , where T D and m x are respectively 
the decoupling temperature and the mass of the WIMPs. 8 ' 10 (Detailed derivations of 
Sv a , and Sv x — including the possible effect of reheating by the e~e+ annihilations — 
can also be found in Ref. 17). Thus, stable and generic caustic formation in the CDM 
infall is inevitable, once density perturbations enter the nonlinear regime. In this 
section, I review the properties of the minimal caustic structure 8-10 that must occur 
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in galactic halos: outer caustics and caustic rings. 
3.1. Outer caustics 

Outer caustics are topological spheres enveloping the galaxies. In catastrophe theory, 
they are known as simple fold (A 2 ) catastrophes. To derive the properties of the 




Fig. 3. Snapshot of the phase space of CDM particles in a galactic halo, r is galactoccntric 
distance and r is radial velocity. The solid line indicates the location of the particles. The dotted 
line corresponds to observer position. Each intersection of the solid and dotted lines corresponds 
to a CDM flow at the observer's location O. Caustics occur wherever the phase space line folds 
back. 

outer caustics, it is instructive to discuss their formation in the phase space point of 
view. 8-10 The small velocity dispersion 5v means that CDM particles lie on a very 
thin 3D sheet in 6D phase space. The thickness of the sheet is the Sv of the particles. 
The flows can be described in terms of the evolution of this sheet. 5 ' 8 Figure 3 depicts 
two dimensional cut (r, r) of the phase space, r being the radial coordinate and r 
being the radial velocity. 8 Before density perturbations enter the nonlinear regime 
of structure formation, there is only one flow (i.e. a single value for velocity) at a 
typical location in physical space. When a large overdensity enters the nonlinear 
regime, the particles in the vicinity of the overdensity fall back onto it. This implies 
that the phase space sheet winds up in phase space wherever overdensity grows. 
The outcome of this process is an odd number of discrete flows at any point in a 
galactic halo; see Fig. 3. Between the first turnaround radius and the location where 
the phase space sheet folds back for the first time — near d the second turnaround 



c "Turnaround" refers to the moments in a particle's history when it has zero radial velocity with 
respect to the galactic center. 

d The ra-th fold is located near the n+l-th turnaround radius; see Fig. 3. If the flow were stationary 
the n-th fold would be located exactly at the n+l-th turnaround radius. 
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radius — there is one flow (associated with particles falling in through the galaxy for 
the first time). Between the first and second folds there are three flows (associated 
with particles falling in through the galaxy for the first and second times, and 
with particles falling out of the galaxy for the first time). Between the second and 
third folds there are five flows and so on. (Wherever the phase space sheet folds 
back, the number of flows in physical space jumps by two naturally). The surface in 
physical space which envelopes the region with the two extra flows, associated with 
the corresponding fold of the phase space sheet, is a CDM caustic. At the fold, the 
phase space sheet is tangent to velocity space, particles pile up there and the physical 
space density diverges as the inverse square root of the distance a to the caustic, 
when it is approached from the side with the two extra flows. This can be shown 
as follows. 17 The 2D phase space mass density near a fold may be represented as 
cPM/dada = MQ(<t)6(<t— (C<t) 2 ), where C and M are positive numbers; see Fig. 3. 
The direction of a coordinate is chosen pointing inward. The density is the integral 
of the phase space density over velocity space 

d{a) = M J d&^^[5(a-ai/C) + 5(& + a^/C)]a-i = AQ{a)a~^ , (1) 

where A is called the fold coefficient. To estimate 16 A, consider the time evolution 
of CDM particles which are falling out of a galactic halo for the nth time and then 
falling back in. Let R n be the nth turn-around radius. We assume that the rotation 
curve of the galaxy is flat near r = R n with time-independent rotation velocity v ro t • 
Then, using the force per unit mass f(r) — —v 2 ot f/r and conservation of energy 
per unit mass E = ^r 2 + V(r), the gravitational potential V(r) is calculated as 
V{r) — v 2 ot In (r/R n ) + E. Thus, the particles have a trajectory r{t) such that 
H = w ro t -\/2 In (R n jr). Here, we neglected the angular momentum of the particles, 
since at the turaround the particles are far from their distance of closest approach 
to the galactic center. Finally, we assume that the flow is stationary. This means 
the number of particles flowing per unit solid angle and per unit time, dN/dildt, 
is independent of t and r, and hence, the caustic is located exactly at the (n + 
l)th turnaround radius R n . The mass density of particles, d n (r), follows from the 
equality: 2dM = d n (r)r 2 dilrdt. The factor of two appears because there are two 
distinct flows, out and in. Near the nth outer caustic where r — R n — a, In (R n /r) ~ 
a/R n . Therefore, d n (r) ~ A n Q{a)a~ 1 / 2 with A n = V2(dM /dndt\ n )v m lRn 3/2 . In 
a self-similar infall, R n and A n are estimated 16 as 

{ft, : n = 1, 2, . . .} =! (240, 120, 90, 70, 60, , , ,) kpc ■ L^'°\ ) (jf) . (2) 

M. : n = 1, 2, . . .} ~ (2, 2, 2, 3, 3, . . .) . , (_Hg_) * (_t) * , („ 

Evidence for the first outer CDM caustic was recently found 12 in the NGC 5846 
group of galaxies. The plot of the surface number density of galaxies was inferred to 
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mark the first outer caustic of radius R\ by an abrupt density drop and a transition 
from a large Sv to a very small 5v, as predicted by the CDM infall models. 22 

3.2. Caustic rings 

Caustic rings were first noted as topological circles in Ref. 8, but their properties 
were left unexplored there. In Rcf. 10 they were studied more precisely as closed 
tubes whose transverse cross-section is a D_4 catastrophe; see Fig. 1. Treatment of 
caustic rings, from the catastrophe theory point of view, is given in Ref. 19. Caustic 
rings are made up of three fold caustics joining at three topological circle of dual 
cusp (A-3) catastrophes. Inside (outside) the caustic rings, there are four (two) 
flows. For simplicity, we consider caustic rings which are axially symmetric about 
the z direction as well as reflection symmetric with respect to z = plane. The 
CDM flow in that case is effectively two dimensional. In galactocentric cylindrical 
coordinates, the flow near such a ring is described by: 

p(Xi,X2) = a + xi 2 -X2 2 , 

z(xi,X2) = 2C(>/p-Xi)X2 , (4) 

where Xi = -\/¥( To — T )' X2 = \f\ a - P 1S distance to the z-axis. The parameters 
t and a label the particles in the flow, t is the time a particle crosses the z = 
plane, a is the declination of the particle, relative to the z = plane, when it was 
at last turnaround, a is the radius of the caustic ring. The constants b, s, u and To 
arc the other characteristics. £ = -^=, and p = \ut^ is the longitudinal dimension 
of the caustic ring cross-section. 6 The number density of particles in physical space 



1 <^ d 2 N _ ^ 1 



2tt P ^ d X id X 2 wvr ' \D 2 (x) 



(5) 

X—Xj(p, z, t) 



where X = (Xi>X2)i an d Xj = (Xi)X2)j ; with j = 1. . .n, are the solutions of 
P = p{XiTX2\t) an d z = z(xi 7 X2' : t)- The number of distinct solutions (flows) n, is 
a function of the location in space p, z, and time t. For the flow near caustic rings, 
the determinant of the Jacobian matrix 

£ 2 = -4C[(xi-VW2) 2 +X2 2 -p/4] . (6) 

Caustic rings occur where D2 vanishes, i.e. where the map (X11X2) ~~ * {Pi z ) is 
singular. Hence, at the caustic, the density diverges in the limit of zero velocity 
dispersion. In the parameter space of X i and X2, the curve for which D 2 = 
is a circle centered at (xi,X2) = (\/p/2,0) with radius ^/p/2. A one-parameter 
representation of this critical circle can be given 19 as: 

Xi = (! + cosV) , X2^^-^^ , (7) 



e If £ = 1, the caustic ring cross-section is a triaxial tricusp. 
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where the angular variable ip e [0, 2tt}. Substituting Eq. (7) into Eqs. (4) yields the 
equations that describe the cross-section of the caustic ring in the (p, z)-plane 

p(ip) = a + - cosi/>(l + cosV'), z(ip) = (- sin?/>(l — costp) . (8) 

Figure 1 shows a plot of this cross-section. The cusps occur at ip = 0, 27r/3,47r/3. 
In a self-similar infall, the caustic ring radii are estimated 9 as 

K : n = 1, 2, . . .} ~ (39, 19.5, 13, 10, 8, . . .)kpc/ #^ ( ( -^-) ,(9) 



0.27 J V h J V220km/s y 

where j max is a parameter proportional to the amount of angular momentum that 
the particles have, h is the present Hubble constant in units of 100 km/(s Mpc) 
and Vrot is the rotation velocity of the galaxy. Because the caustic rings lie close 
to the galactic plane, they cause bumps in the rotation curve, at the locations of 
the rings. In a study of 32 extended and well-measured external galactic rotation 
curves, evidence was found for the law given in Eq. (9). 13 In the rotation curve of 
the Milky Way, the locations of eight sharp rises fit the prediction of the self-similar 
model at the 3% level. 14 Moreover, a triangular feature in the Infrared Astronomy 
Satellite (IRAS) map of the Milky Way plane was interpreted as the imprint on 
baryonic matter of the fifth caustic ring. 14 This nearby ring is 55 pc away from us. 
The feature is correctly oriented with respect to the galactic plane and the galactic 
center. Its location is consistent with the bump between 8.28 kpc and 8.43 kpc in 
the rotation curve. The probability that the coincidence in position of the feature 
with a rise in the rotation curve is less than 10~ 3 . 

The density profile near the caustic ring along the direction a perpendicular to 
the surface is d(tp,cr) = A(i(j)0(a)a^ 1 ' >2 , where the fold coefficient 



., ., d 2 M 2( cosa(V0 nn , 



with 19 



(1 + 2cosi/i)tan^ 



(1 +cosV>) 2 +C 2 sin 2 tp . (11) 



The fold coefficients are estimated as a function of location on the ring as 
{A n (V>):n=l,2,. . .}~(5,6,6,7,8,...)-l(r 



cm 2 kpc 2 



where 19 



' 1 °1) (o't) (220km s> (]2) 



■FnWO - , , — = • (13) 

yf\ (1 + 2 cos V) tan f | [(1 + cos V) 2 + CI sin 2 V] * 
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Due to the symmetry of the tricusp with respect to the z — axis, we have ^ r n ('i?) = 
Fn{— $) an d T n {^f ± -d) = f n (y ^ iJ). For the typical values p n = 0.1 a n and 
( n = 1, near the outer cusp ^"(i?) — 0.525 ?9~2. Near the nonplanar cusps T{^- ± 

11 

1?) ~ 0.585 ?9 _ 2. This divergence of the fold coefficients at the cusps is cut off 
when Sv ^ 0, because, the location of the caustic surface gets smeared over some 
distance Sx and the cusps are smoothed out. For CDM caustics in galactic halos, 
Sx and Sv are related 9 ' 15 ' 16 as Sx ~ R 5v/v where v is the order of magnitude 
of the velocity of the particles in the flow and R is the distance scale over which 
that flow turns around. For a galaxy like our own, v = 500 km/s and R = 200 
kpc are typical orders of magnitude. As a result of primordial velocity dispersion, 
axion and WIMP caustics in galactic halos are typically smeared over Sx a ~ 6 • 
10 4 km (lO~ 5 eV/m a ) 5 ^ 6 and 8x x <~ 10 10 km (GeV/m x ) 1 ^ 2 , respectively. However, a 
CDM flow may have an effective velocity dispersion Sx c s which is larger than its 
primordial velocity dispersion. Effective velocity dispersion occurs when the sheet 
on which the dark matter particles lie in phase space is wrapped up on sub-scales 
that are small compared to the galaxy as a whole. Little is known about the size of 
the Sx e g in galactic halos. However, the widths of the bumps at the ring locations in 
the rotation curve of the Milky Way and the sharpness of the edges of the triangular 
feature in the IRAS map were used to deduce an upper bound of 15-20 pc on Sx e s 
at most. 14 

I use the primordial smearing distances Sx a and 5x x to estimate the upper limits 
for the fold coefficient and for the lensing effects at the cusps of caustic rings as 
follows. I set the Sv of the particle species equal to zero, and then calculate the 
fold coefficients and the lensing effects (Sect. 4) on the caustic surface at a point 
whose transverse distances Ap and Az to the nearby cusp — that would occur in the 
limit Sv = — are slightly greater than the primordial smearing out distance Sx a or 
Sx x . These points correspond to the locations that are respectively Aip = and 
Aip = j^r radian away from the cusps in the axionic and WIMPic cases, respectively. 
Similarly, I use the upper bound for the smearing distance Sx e s to estimate the lower 
limits: I choose a point which is about the maximum smearing distance Sx c s away 
from the nearby cusp and calculate the fold coefficients and the lensing effects there. 
Such points differ by Aip = radians from the cusp locations. The calculations 19 
show that, for the axion caustic rings, the factor ^(tp) may change the estimates 
given for the A n in Eq. (12) between 0.8 and 26 times at the outer cusp, and between 
0.9 and 29 times at the nonplanar cusps. The upper bounds of A n for the WIMP 
caustic rings, on the other hand, are ten times less then the upper bounds for the 
axionic rings. 

4. Gravitational Lensing by the CDM Caustics 

Matter curves the spacetime. As a result, light rays coming from a background 
source are deflected due to the presence of a forefront object. This effect is known 
as gravitational lensing. Choosing the y-axis in the direction of propagation of light, 
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the angular shift in the position of the source £s = (Jisxi ^Sz) due to a lens is given 
by A£ = £j — £s = Vg 7 $(£j), where = (£j x , £/ z ) is the position of the image with 
lens present. The 2D potential $ is solved from the Poisson equation V| 7 $ = 2£/£ c , 
where column density £ and the critical surface density E c are defined as: 

= J dyd{DU Ix ,y,DU Iz ) , 

c2D S no.-, 8 f A? 



— 4^ Mm -°- 347 c^^ G P C J- <"> 

Here, £>x, and D$ are the distances of the observer to the lens and to the source, 
respectively. Dls is the distance of the source to the lens. The image distortion and 
magnification are given by the Jacobian matrix of the map £s(£j)- Because grav- 
itational lensing does not change surface brightness, the magnification M. is the 
ratio of image area to source area. Hence, A4 = 1/| det {d£si/d£ij)\, which yields, 
to first order, M. = 1 + 2£/£ c . Gravitational lensing has proven useful in revealing 
the massive compact halo objects (MACHOs) in galaxies and constraining the mass 
distribution in galaxy clusters. For some recent interesting lensing applications see 
Refs. 23. Dark matter caustics have calculable lensing signatures. 15-19 The caustic 
ring model 8-10 precisely predicts the density and the geometry of the CDM distri- 
bution in caustic neighborhood. In Rcf. 16, we derived the lensing equations for the 
line of sights that are parallel to the galactic plane of a caustic ring and near tangent 
to its surface in some specific cases. In Ref. 19, I obtained the lensing equations for 
the same line of sights, but near tangent to a caustic ring anywhere on the surface 
and estimated the lensing effects near the cusps of axion and WIMP caustic rings. 
Here, I summarize these effects. 

The gravitational lensing effects of a caustic are largest when the line of sight 
is near tangent to the surface, because the contrast in column density X is largest 
there. For such line of sights, the associated fold of outer caustics, in general, and 
of caustic rings, in the regions where < ip < 2n/3 and 4n/3 < ip < 2ir, are curved 
towards the side with the two extra flows; see Fig. 4. f In these cases, the image shift 
and magnification are given 16 as 

A£ = &-& = »?£/e(-&), M^§^ = i + v e(-Ci) + o( v 2 ) 1 (is) 



where 



2nA^2\R\ 

V = ^ ■ (16) 



Here R is the curvature radius of the surface along the direction associated with 
the chosen line of sight. Hence, when the line of sight of a moving source crosses 



f If the line of sight is near tangent to the ring surface at ip € the associated fold is curved 

away from the side with the two extra flows. The expected effects are not as large in that case. 16 ' 19 



5, 2008 15:11 WSPC/INSTRUCTION FILE Oncmli-cusplens 



12 V. K. Onemli 




B 




Fig. 4. Lcnsing by a caustic ring for a line of sight near tangent to the surface at a given point tp* 
close to the outer cusp. The line of sight lies in the z = z(?/>*)-plane. A) Side view in the direction 
of the line of sight. The latter is represented by the dot near x = z = 0. B) Top view. 



the caustic, the component of its apparent velocity perpendicular to the surface and 
the magnification of the image changes abruptly. Both affects are of order rj. 

The r\ n for the outer caustics are estimated, 16 using Eqs. (2) and (3) in Eq. (16), 

as 

{ Vn :n = l,2... .}~(7, 6, 6, 5, 5, . . .) • lO" 3 

X (d s Gpc) (220km/s) (oj) ' 

Lensing effects are always maximum when the lens is situated half-way between the 
source and the observer. Also, D$ should be as large as possible to get the largest 
effects. Therefore, in my estimates, I will assume that the source is at cosmological 
distances, e.g. 2D L = 2D LS = D s = 1 Gpc. For the outer caustics, this yields 
a magnification of order 10~ 3 which can not be detected by present instruments. 
However, the images of extended sources may be modified in recognizable ways. 15,16 
In particular, if a straight jet makes angle a with the normal, it appears bent by 
an angle ^sin(2a)/2 where its line of sight crosses a fold of an outer caustic. 16 
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For the caustic rings, A is given in Eq. (10). Along the chosen line of sights, the 
curvature radius of the associated fold 19 

J R = -C _1 Sign(cos^-cos2^)[a+|cos^(l + cos-0)](C 2 +cot 2 ^)3 . (18) 
Therefore 19 

WMcc^ 
S c dudr b^/ap 

where 



GW - < / ^ r • (20) 

y [1 + £ cos V(l + cos V)] I (1 + 2 cos V) tan f | 
The JfoWO for the rings are estimated 19 as 

{?7„(V) : n - 1,2,...} ~ (7,6,6,5,5,...) • lO" 2 ^) ° L ^ 



D s Gpc 
0.27 \/ h\f v rot 



Jma X A0-7A 22 0km/s 



(21) 



In the limit of zero velocity dispersion Q(ip) diverges at the cusps. G("&) — 0.8 # _1 , 
whereas (?(27r/3 — t?) ~ 0.6 i9~ 1//2 . Hence, the effects are large near the cusps, in 
particular near the outer cusp. For the line of sights near tangent to a cosmological 
ring surface in the region — 7r/27 < tp < n/27, the magnification is more than 10%. 
Effects of that order can be observed by present instruments. 24 

The lensing effects of the cusps of axion and WIMP caustic rings may be 
estimated 19 by turning the 5v of the particle species on, as described in Sect. 
3.2. At the outer cusp of a cosmological axion caustic ring, the range of r\ is con- 
strained between 7j(±y^) ~ 2.9 • 10 -2 and rj(±j^) ~ 27.9. This implies a magnifi- 
cation between 3% and 2800%. At the outer cusp of a cosmological WIMP caustic 
ring, on the other hand, the range is constrained between ^(i^) ~ 0.029 and 
i](±j^) ~ 0.28. This implies a magnification between 3% and 28% at the outer 
cusp. The lower bound of r] at the nonplanar cusps of cosmological caustic rings 
is rj(^- — ~ 0.016. For a cosmological axion caustic ring, the upper bound is 
r}(^§- — 7550) — 0.46. Thus, the magnification can be between 2% and 46% at the 
nonplanar cusps. (Twenty five MACHOs were discovered 24 by magnifications that 
range between 12% and 46%). For a cosmological WIMP caustic ring, on the other 
hand, — jf) — 0.046, hence the magnification can be between 2% and 5% at 
the nonplanar cusps. 

Pointlike background sources may probe the CDM caustics, as they cross behind 
the cusps of a foreground halo. The time scale St for the brightness to change is 
about the transit time of a background source across the thickness (i.e. smearing out 
distance) Sx of the caustic edges. For axion caustics St a is about an hour, whereas 
for WIMP caustics 8t x is about a year. 15 A strategy to detect the caustics may 
be a MACHO style experiment: monitoring large number of pointlike background 
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sources behind the disks of foreground halos and watching for their edge crossings 
near the cusps. The number of outer caustics across the face of a halo is about the 
number of orbits (discrete flows) since it is formed, 5,15 say N s=s 10 to 100. The 
number of caustic rings are of the same order, because, the very same flows that 
generate the outer caustics, generate the caustic rings. 8 Caustic ring radii a n ~ 40/n 
kpc. Each background source traverses such a distance scale about 10 8 /n years. So 
monitoring 10 8 stars for a year will show about 20 to 200 crossings of the ring 
edges. Only 6% of the ring surface (near the cusps) can magnify the sources more 
than 10%, hence about 1 to 10 cusp crossings may be detected in a year. Note 
that I estimated the prospects considering the minimal cusp structure predicted by 
the caustic ring model. Cusps are generic and common throughout the halos. More 
background sources may easily pass behind the generic cusps of a foreground halo. 

5. Conclusions 

Caustics are common features in the propagation of light, because photons are 
collisionless and flow of light from a point source has zero velocity dispersion. In 
the occasions where these two conditions are satisfied, caustics in the flow of the 
ordinary luminous matter are also observed. Cold dark matter (CDM) flow onto 
isolated galaxies has both of the properties; the particles interact only weakly and 
the flow has a tiny velocity dispersion. Therefore, caustic formation in the halos 
of CDM is inevitable. Gravitational lensing effects of the cusps at cosmological 
distances may be detected by present instruments. Caustic ring model of galactic 
halos predicts the geometry and density near the cusps of the rings. I derived the 
lensing equations for the line of sights near tangent to a caustic ring anywhere 
on the surface, and estimated the lensing effects near its cusps. For a cosmological 
axion caustic ring, the magnification may range between 3% and 2800% at the outer 
cusp, and between 2% and 46% at the nonplanar cusps. For a cosmological WIMP 
caustic ring, on the other hand, the magnification may range between 3% and 28% 
at the outer cusp, and between 2% and 5% at the nonplanar cusps. As pointlikc 
background sources pass behind the caustics of foreground halos, the time scales 
for brightness to change is about an hour for the axion folds and about a year for 
the WIMP folds. Hence, depending on the strength of the observed effect and the 
time scale for brightness change, it may even be possible to discriminate between 
the axions and the WIMPs. 
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g The sharp rises in the rotation curve of the Milky Way imply the existence of 13 caustic rings. 
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